Inelastic light scattering from a Mott insulator 
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We propose to use Bragg spectroscopy to measure the excitation spectrum of the Mott insulator 
state of an atomic Bose gas in an optical lattice. We calculate the structure factor of the Mott 
insulator taking into account both the selfenergy corrections of the atoms and the corresponding 
dressing of the atom-photon interaction. We determine the scattering rate of photons in the stim- 
ulated Raman transition and show that by measuring this scattering rate in an experiment, in 
particular the excitation gap of the Mott insulator can be determined. 

PACS numbers: 03.75.Hh, 67.40.-w, 32.80.Pj, 39.25+k 



Introduction. — A Bose-Einstein condensate in an op- 
tical lattice is a powerful new tool to investigate strongly 
correlated Bose gases El @. In particular, the experi- 
ment by Greiner et al. [3j has shown that it is possible 
to achieve a quantum phase transition from a superfluid 
to a Mott-insulating phase by starting from an atomic 
Bose-Einstein condensate and loading this system into 
an optical lattice. This phase transition was predicted 
to occur in the Bose-Hubbard model by Fisher et al. 0, 
and Jaksch et al. [f| were the first to make the crucial ob- 
servation that the Bose-Hubbard model can be applied 
to bosonic atoms in an optical lattice. A mean-field the- 
ory that describes the two phases of the Bose-Hubbard 
model was developed by van Oosten et al. [||. 

An important advantage of using atoms in an opti- 
cal lattice to study the Bose-Hubbard model, is that the 
system is free from disorder, which makes it possible to 
make very accurate predictions and measurements. A 
good example of such a high-precision measurement is 
Bragg spectroscopy. This technique has already been 
used to coherently split a Bose-Einstein condensate into 
two momentum components 0, to measure the excita- 
tion spectrum of a trapped Bose-Einstein condensate || , 
and to measure the light-shifted energy levels of an atom 
in an optical lattice Q. Here we propose to use Bragg 
spectroscopy to measure the excitation spectrum of the 
Mott-insulator state. In particular, one can in this way 
determine the value of the particle-hole gap in the excita- 
tion spectrum and study the behaviour of this gap as the 
system approaches the quantum critical point. Note that 
the excitation spectrum as obtained using Bragg spec- 
troscopy will not yield what is generally referred to as 
the Mott gap, because this gap is associated with single- 
particle excitations. The value of the particle-hole gap is 
an interesting quantity in the study of quantum critical 
phenomena, but it is also very important for the practi- 
cal application of these systems in quantum information 
processing, since the gap determines the fidelity of the 
Mott state. It is important to realize that the system of 
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FIG. 1: (a) Setup for the proposed experiment, (b) Parti- 
cle and hole dispersions in units of the tunneling parameter 
in a one dimensional lattice, for U/zt — 6. The horizontal 
arrow indicates absorption of momentum, the vertical arrow 
absorption of energy. 



a Bose-Einstein condensate in an optical lattice is more 
complicated than the above mentioned systems, because 
in this case many-body effects and strong correlations 
have to be taken into account. In a Bragg spectroscopy 
experiment, two laser beams are used to make excitations 
in the system, as shown in Fig. QJa) . The two lasers both 
have a large detuning with respect to an optical transition 
in the atoms so that spontaneous emission is suppressed. 
However, the relative detuning can be very small. When 
an atom absorbs a photon from beam two and is stim- 
ulated to emit a photon into beam one, the atom un- 
dergoes a change of momentum hq = Kk.2 — Kki and a 
change of energy tku = T10J2 — Tiuj\. In principle any op- 
tical transition could be used, but here we use the same 
transition that is employed to create the lattice poten- 
tial. This means that the magnitude of the momentum 
is given by Tiq = 2hk p \ l sin(0/2), where to a good approx- 
imation ftfcph = 2tt%/\ is the photon momentum of both 
the lasers, A is equal to the wavelength of the lattice laser 
light and is the angle between the two laser beams. By 
varying the angle between the two laser beams, any mo- 
mentum between zero and 2hk p h can be transferred and 
by varying the relative detuning between the beams, the 
amount of energy that is transferred to the system can 
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be controlled. 

Calculating the scattering rate for a given momentum 
hq and energy Tiuj roughly speaking involves counting the 
number of ways in which the requirements of momentum 
and energy conservation can be met. To illustrate this 
process, we draw in Fig.QJb) the quasiparticle and quasi- 
hole dispersions in the Mott insulator as is common 
in solid-state physics. The horizontal and vertical arrows 
in the figure indicate the transfer of momentum and en- 
ergy respectively. Since energy is deposited in the system, 
this scattering rate can be measured in a traploss exper- 
iment, or by determining the increase in temperature of 
the atoms. 

The scattering rate. — To calculate the desired two- 
photon scattering rate we use Fermi's Golden Rule. In 
linear response, this can be expressed as /(q, u>) — 
— 21m [n(q, u>)] /h, where n(q, uj) is the polarizibility 
of the medium. The polarizibility can be written as 
II(q, to) — (hfl/2) 2 x(q, w), with fl the effective Rabi fre- 
quency for the two-photon process and x the susceptibil- 
ity. The retarded susceptibility is given by 

X (q,w) = - jjdx p dt'e-^-^ 

x ([at(x ) t / )a(x,t'),a + (O,0)a(O,0)]) ) (1) 

with V the volume and cr (x, if) and d(x, t') creation and 
annihilation operators of the atoms. Because the atoms 
are in an optical lattice, we can expand the field operators 
in terms of the Wannier states of the lattice, which yields 
an expression in terms of creation and annihilation oper- 
ators for every lattice site. As mentioned previously, the 
Hamiltonian of the system then equals the Bose-Hubbard 
Hamiltonian with a tunneling amplitude t, an on-site in- 
teraction energy U and a chemical potential /z. Using the 
decoupling approach described in Ref. Q, we can write 
the atomic propagator in the Mott-insulator phase as 



Z(k) 



1 - Z(k) 



-hcu+ + e ( iP(k) ~hu+ + e qh (k) 



(2) 



where the probabilities Z(k) and 1 — Z(k) account for 
the fact that an atomic excitation contains both quasi- 
particle and quasihole contributions. The notation fiu> + 
is shorthand for hco + i£ with £ J. 0. The dispersions for 
the quasiparticle and quasihole excitations are given by: 

e^ h (k) = - M +|(2iVo-l) + i(e(k)±Mk))(3) 

where A^o is the filling fraction of the lattice and the 
function e(k) = — t Ylj=i cos 27rfcj corresponds to the lat- 
tice dispersion in the experimentally relevant case of a 
regular square lattice. The momentum hit is here and 
from now on always written in units of 2hk p ^, which 
means that the first Brillioun zone runs from kj — 
— 1/2 to 1/2. The energy ?kj(k) is given by KujQl) = 
y/U 2 + (4N + 2)J7e(k) + e(k) 2 and the probability Z(k) 



is given by Z(k) = (U(2N + 1) + e(k) + fiw(k))/2/iw(k). 
Using the Greens function in Eq. (|2J), we find in first ap- 
proximation x (q,w) = t(q) (x<{ (q,w) -X+(-q,-w)), 
where t(q) is a geometric factor that involves the appro- 
priate overlap integral of the relevant Wannier functions 
(this will be discussed in Ref. 0). Denoting integration 
over the first Brillouin zone as / 1BZ ,the contribution due 
to the creation of a particle-hole pair is given by 

0/ a 1 / P(k,k + q,^) 



e q P(k- 



e qh (k)' 



(4) 

and the time-reverse process can be written as 
X^(q> = X+( — This equation contains the 
probability P(k, k + q, u) = (1 - Z(k)) Z(k + q) for the 
creation of a hole with momentum k and a particle with 
momentum k + q, and an energy denominator that is as- 
sociated with the energy cost e qp (k + q) — e qh (k) of that 
process. This can readily be verified by taking the imag- 
inary part of the susceptibility, which is proportional to 
J mz dkP(k, k + q, uj)S{huj - e qp (k + q) + e qh (k)) and can 
be understood as Fermi's Golden Rule. The actual com- 
putation of the above integral is too complicated to do 
analytically, so that we have to resort to numerical meth- 
ods. We achieve this by calculating the imaginary part 
of Eq. , which roughly corresponds to integrating over 
the surface in the Brillouin zone where the energy denom- 
inator vanishes. In practice, this amounts to numerically 
finding the poles of the expression and determining their 
residue. The real part is calculated from the imaginary 
part using a Kramers-Kronig relation. 

However, the results that one would obtain in this 
manner do not obey particle conservation. Physically, 
a Raman process with momentum q couples to a den- 
sity fluctuation p(q). For zero- momentum transfer, p(0) 
corresponds to the total number of particles and fluctua- 
tions are impossible due to particle-number conservation. 
If we compute the imaginary part of Eq. (j^J for q = we 
find a spectrum which is nonzero, which means that this 
approach is not sufficiently accurate. The problem is due 




FIG. 2: Schematic representation of (a) Eq. (gj and (b) 
Eq. © 

to the fact that in Eq. J5J not the bare atomic propagator 
is used, but a dressed propagator which contains a large 
self-energy correction given by 

N (N + 1)U 2 



aS(k,cj) = 2N U 



(5) 



fkj + U + fi 

The first term on the right-hand side is the Hartree-Fock 
contribution, which is also present in a Bose-Einstein con- 
densate. The second contribution is due to the correla- 
tions in the Mott insulator. Essentially this means that 
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an atom moving through the Mott insulating background 
is dressed by all the other atoms. As is known from 
quantum field theory one has to be careful when 
applying self-energy corrections to the calculation of the 
susceptibility, because in general these corrections do not 
obey the required conservation laws (in this case particle- 
number conservation). Using field-theoretical methods, 
we can derive so-called Ward identities, that show that 
every self-energy correction requires a corresponding ver- 
tex correction in order to restore the conservation laws. 
Physically this means, that if the atom is dressed, we also 
have to dress the atom-photon coupling. Diagrammati- 
cally this is illustrated in Fig. This situation is very 
analogous to the situation in a superconductor, where the 
naive BCS calculation of the electro-magnetic response 
is not g aug e invariant and a more involved approach is 
needed 0|- Using the relevant Ward identity we 
can derive that the intuitive probability function given 
above, has to be replaced by 



P(k,k + q,^) = 



2fto;-eqP(k + q) + e qh (k) 

?ia;(k + q) + Tiw(k) 
(Z(k + q)-Z(k)). 



(0) 



Note that the probability now vanishes when q — ► 0, so 
that particle conservation is indeed no longer violated. 
In fact, we can show that for small q and Hlo just above 
threshold P oc g 2 /Ag, where Ao is the gap for particle- 
hole excitations. 
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FIG. 3: Real and imaginary parts of the susceptibility for 
U/zt = 10 and q = 0.10,0.14,0.18,0.20 along a lattice di- 
rection, in two dimensions. The dotted line in the bottom 
figure is the result for q = 0.001 multiplied by 250 to show 
the behaviour for small q. 

Results. — In Figs. and 0] the result of a numer- 
ical integration is shown in two and three dimensions, 
respectively. Both calculations have been carried out for 
a regular square lattice and the momentum q is chosen 
a principal lattice direction. All energies in the following 
figures are given in units zt, where z is the coordination 
number of the lattice. 



The imaginary part of Fig. [3] clearly shows singulari- 
ties around %u> = U . These singularities are due to the 
fact that there are saddle points in the dispersion and 
that a saddle point in the dispersion causes an integrable 
singularity in the density of states. These are so-called 
van Hove singularities [l3|. It is interesting to see, that 
the van Hove singularities split up as the momentum is 
increased, which is caused by the fact that the saddle- 
point energy in the direction of q and the saddle-point 
energy in the orthogonal direction(s) are shifted by dif- 
ferent amounts. This is also visible in Fig.^J However, it 
is less clear in this case, because the van Hove singulari- 
ties are more smeared out in three dimensions. Also, the 
opening of the threshold for the two-photon absorption 
in the three dimensional case is far less steep than in the 
two-dimensional case. To investigate possible collective 




FIG. 4: Real and imaginary parts of the susceptibility for 
U/zt = 10 and q = 0.10,0.14,0.18,0.20 along a lattice di- 
rection, in three dimensions. The dotted line in the bottom 
figure is the result for q = 0.001 multiplied by 250 to show 
the behaviour for small q. 

modes in this system, we determined higher-order cor- 
rections in the random-phase approximation (RPA). It 
can be shown that in RPA the susceptibility is given by 
X(q,w) = x°(q,w)/(l - Ux°(<1,uj))). This means that 
there is a resonance in the scattering rate when the real 
part of X°(q, oj) is equal to 1/U. However, as can be 
seen from Figs. |21 and 01 the real parts in both cases are 
rather small compared to 1/U and in practice, including 
the RPA denominator does not qualitatively change our 
previous results. 

In Fig. we plot the imaginary part of xot for a range 
of values for the coupling constant U/zt, and for a fixed 
momentum q = 0.10. We see, that the threshold be- 
haviour becomes steeper as we approach the critical value 
of U c / zt as 5.83. We also see that there remains a nonzero 
gap when U = U c . This is due to the fact that we are 
not considering a zero momentum excitation, due to the 
reasons given above. In the inset of Fig. we plot this 
gap A q as a function of U/zt. For large U the gap grows 
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linearly with U and for U close to U c , the gap closes more 
rapidly. In the case of q = the gap would in our mean- 
field approximation close as \JU — U c when U I U c , but 
for small nonzero q it closes as y/U — U c + rjq 4 , where 
the factor 77 is a positive function of U c and t. 
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FIG. 5: Imaginary part of the susceptibility in a three di- 
mensional lattice for q = 0.10 along a lattice direction and 
U/zt = 5.83,7,8,9, 10 

Discussion. — In summary, we have proposed a means 
of studying Mott insulators in optical lattices, using the 
relatively well-known technique of Bragg spectroscopy. 
We have presented spectra that can be measured directly 
by trap loss or heating measurements. In a recent experi- 
ment by Stoferle et al. Q the authors use a setup where 
the laser beams are perfectly counterpropagating, which 
corresponds to a quasi-momentum transfer of zero. As 
we have argued above, there should be no scattering in 
that case and the signal can only be due to nonlinear re- 
sponse or to the fact that the system is inhomogeneous 
and of finite size. We have found that by measuring the 
threshold behaviour of the two-photon scattering rate at 
various quasi-momenta, it is possible to determine the 
gap by extrapolation. We have shown that for a theoret- 



ical description of Bragg spectroscopy on the Mott insu- 
lator it is absolutely essential to dress the photon-atom 
coupling, which is in a way unexpected, as the corrections 
are zero in the case of an harmonically trapped gas. As 
a result it turns out that although it is common to use 
the language of solid-state physics to describe these sys- 
tems, the physics is qualitatively very different due to the 
many-body effects. 
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